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In this paper the constitutive relations for moving media with homogeneous 
and isotropic electric and magnetic properties are presented as the connections 
between the generalized magnetization-polarization bivector M. and the elec- 
tromagnetic field F. Using the decompositions of F and Ai, it is shown how 
the polarization vector P(x) and the magnetization vector M(x) depend on 
E, B and two different velocity vectors, u - the bulk velocity vector of the 
medium, and v - the velocity vector of the observers who measure E and B 
fields. These constitutive relations with four-dimensional geometric quantities, 
which correctly transform under the Lorentz transformations (LT), are com- 
pared with Minkowski's constitutive relations with the 3-vectors and several 
essential differences are pointed out. They are caused by the fact that, contrary 
to the general opinion, the usual transformations of the 3-vectors E, B, P, M, 
etc. are not the LT. The physical explanation is presented for the existence 
of the magnetoelectric effect in moving media that essentially differs from the 
traditional one. 
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1. Introduction 

In this paper the constitutive relations for moving media with homogeneous and 
isotropic electric and magnetic properties are presented using the abstract four- 
dimensional (4D) geometric quantities and their representations in the standard 
basis. These results are compared with Minkowski's constitutive relations 1 with 
the 3-vectors and with Minkowski's constitutive relations that are obtained by 
means of exterior forms. 2 The paper is organized as follows. 

First, in Sec. 2, we present a review of some previous results that are impor- 
tant for the theory presented here. In the recent paper, 3 a formulation of the field 
equations for moving media is developed by the generalization of an axiomatic 
geometric formulation of the electromagnetism in vacuum. 4 As mentioned in 
Ref. 3 almost the entire physics literature deals with the electromagnetic exci- 
tation tensor H, i.e. with the electric D and magnetic H excitations (see Eq. 
([5])) and the constitutive relations refer to the connections between them and 
F, i.e. E and B, respectively. But, as shown in Ref. 3, it is physically better 
founded to formulate the field equations for moving media in terms of F and 
the generalized magnetization-polarization bivector M instead of, as usual, F 



1 



and H. In Sec. 2, the decompositions of F (JJ and M are presented. F 
is decomposed into vectors E, B and v - the velocity vector of the observers 
who measure E and B fields. M is decomposed into the polarization vector P, 
the magnetization vector M and u - the bulk velocity vector of the medium. In 
Ref. 3, the field equations are written in terms of F and M. and also in terms 
of vectors E, B, P, M and the velocity vectors u and v. These field equations 
are also quoted in Sec. 2. Furthermore, in Sec. 2, the usual transformations 
(UT) of the 3-vectors E and B, ©, and of P and M, ([7]), are written and 
their difference relative to the Lorentz transformations (LT) of vectors, as 4D 
geometric quantities, e.g., the electric field vector E, ©, is pointed out. 

In Sec. 3, we formulate the constitutive relations as the relations between 
M and F, Eqs. (fTTj) and (fl2j) . Then, using the decompositions of F fll} and M 
(J3]> we get from ((TTJ) and ([12]) how P(x), (fT3]), and M(x), dHJ), depend on E, 
B and two different velocity vectors, u and v. The constitutive relations (|13|) 
and (TTJI are the basic relations that are obtained in this paper and they are not 
reported in previous approaches. The last term in (fT3"|) and (fT4")) describes the 
magnetoelectric effect in a moving dielectric in a new way. 

In Sec. 4, we have represented all 4D geometric quantities from (|T3| and (fT4| 
in the standard basis in order to compare them with some usual formulations. 
This procedure yields Eqs. (ITT)) and (TTS)) . 

In Sec. 5, Minkowski's constitutive relations with the 3-vectors (f2"2"| are 
quoted. The equations (|22p are considered to be the fundamental constitutive 
equations for moving media in the whole physics community. In Sec. 5.1, the 
relations (|22|) are written in equivalent forms as constitutive equations which 
explicitly express the 3-vectors P and M as the functions of the 3-vectors E and 
B, (|23]) and (f24|) . or ([25]) . These forms of Minkowski's constitutive relations are 
compared with our relations (TT7)) and (ITS)) , i.e., with Eqs. (flU)) and (|2"Tj) . and 
several essential differences are pointed out. It is argued that these differences 
appear since Minkowski's constitutive relations are with the 3-vectors and they 
are derived using the UT and not the LT. In Sec. 5.2, it is shown that the 
same differences remain in the low velocity limit. In Sec. 5.3, it is presented the 
comparison with Minkowski's constitutive relations that are obtained by means 
of exterior forms. 2 It is shown that the constitutive relations with exterior forms 
from Ref. 2 are completely equivalent to Minkowski's constitutive relations with 
the 3-vectors and, accordingly, they also differ from the relations obtained in 
this paper. 

In Sec. 6, Eq. (|28]) represents the interaction term in the Lagrangian for the 
interaction between the electromagnetic field F and the dipole moment bivector 
D, whereas Eq. (|29l) is its low velocity limit. The last two terms in (l28l) . or ([29]) , 
contain the direct interaction of E with the magnetic dipole moment vector m 
and B with the electric dipole moment vector d. These terms give the physical 
explanation for the existence of the magnetoelectric effect in moving media. 
That explanation markedly differs from the traditional one. 

In Sec. 7, some remarks are given, which refer to the general constitutive 
relations. 

In Sec. 8, the conclusions are presented. 
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2. A Short Review of Some Previous Results 

We shall deal with 4D geometric quantities, i.e. in the geometric algebra for- 
malism. For the exposition of the geometric algebra see Ref. 5. The genera- 
tors of the spacetime algebra are four basis vectors {7^} , fj. — 0...3, satisfying 

= T\ia> = diag(^ ). This basis, the standard basis, is a right-handed 

orthonormal frame of vectors in the Minkowski spacetime M 4 with 70 in the 
forward light cone, 7q = 1 and 7^ = — 1 (A; = 1,2,3). The standard basis {7^} 
corresponds to Einstein's system of coordinates in which the Einstein synchro- 
nization of distant clocks 6 and Cartesian space coordinates x l are used in the 
chosen inertial frame of reference. 

First, we briefly expose the formulation of the field equations from Ref. 

3. As shown in Ref. 3, the field equations (5), d(e F + M) = j {c) /c; d ■ 
(eoF + M) = j( c ^/c, d A F — 0, can be taken as the primary equations for 
the electromagnetism in moving media. The bivector F = F(x) represents the 
electromagnetic field and, as shown in Ref. 4, it can be taken as the primary 
quantity for the whole electromagnetism. is the conduction current density 
of the free charges. M. is the generalized magnetization-polarization bivector 
M = M(x), which is connected with the magnetization-polarization current 
density of the bound charges j^ M ^ = —cdM = —cd ■ M.. (According to Eq. 
(3), 3 the total current density vector j can be decomposed as j — + j^ M ^ ■) 
The field equation with sources is written in the 'source representation' in Eq. 
(7) 3 d ■ EqF = jc — d ■ M.; the sources of F are the true currents j( c ' and 
the magnetization-polarization current density d ■ A4. In most materials M. is 
a function of the field F and this dependence is determined by the constitutive 
relations. In all previous formulations of the electromagnetism in media (at 
rest, or moving), the electromagnetic excitation tensor is introduced as H = 
SqF + A4. Then, the constitutive relations refer to the connections between H 
and F. However, as discussed in Ref. 3, physically different kind of entities 
are mixed in that definition for H; an electromagnetic field F and a matter 
field Ai. Furthermore, in general, two different velocity vectors, v - the velocity 
of the observers and u - the velocity of the moving medium, enter into the 
decompositions of F and M, respectively. For F, that decomposition is given 
byEq. (10) 3 

F = EAv/c+(IcB)-v/c, (1) 

where the electric and magnetic fields are represented by vectors E(x) and B(x) 
and / is the unit pseudoscalar. There is no rest frame for the field F, that is, 
for E and B, and therefore the vector v in the decomposition ([1]) is interpreted 
as the velocity vector of the observers who measure E and B fields. Then E(x) 
and B(x) are defined with respect to v, i.e., with respect to the observer, as 

E = F-v/c, B = —(l/c)I(F Av/c). (2) 

It also holds that E-v = B -v = 0; both E and B are space-like vectors and they 
depend not only on F but on v as well. Similarly, in Eq. (12), 3 the bivector 
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M(x) is decomposed into two vectors, the polarization vector P(x) and the 
magnetization vector M(x) and the unit time-like vector u/c 



M = P A u/c + (MI) ■ u/c 2 . 



(3) 



There is the rest frame for a medium, i.e., for M., or P and M, and therefore 
the vector u in the decomposition ([3]) is identified with bulk velocity vector of 
the medium in spacetime. Then, P(x) and M(x) are defined with respect to u 
as 



and it holds that P ■ u = M ■ u = 0. As in the case with F, it is visible from (0| 
that P and M depend not only on A4 but on u as well. 

Usually, only the velocity vector u of the moving medium is taken into ac- 
count, or the case u = v is considered, i.e., it is supposed that the observer 
frame is comoving with medium, or both decompositions ([1]) and ([3]) are made 
with the same velocity vector, either u or v. s Such assumptions enable the in- 
troduction of the electromagnetic excitation bivector and, by using ([T]) and 
([3]), one finds the decomposition of % into the electric and magnetic excitations 
(other names of which are 'electric displacement' and 'magnetic field intensity') 



(Eq. (14) 3 ), where, as usual, the electric displacement vector D = eqE + P and 
the magnetic field intensity vector H = (1/ pio)B — M are introduced. 

In Ref. 3, inserting the decompositions of F and A4 (Eqs. ((T|) and (J3J here) 
into the field equation, the general form of the field equation for a magnetized 
and polarized moving medium with E(x), B(x), P(x) and M(x) (the Amperian 
form) is obtained in the form of Eq. (15), i.e., (16) (the vector part, with 
sources), or (18) in the 'source representation' d ■ {£q[E A v/c + (IB) ■ v]} — 
jW/c - d ■ [P A u/c + (l/c 2 )(MI) ■ u], and (17) (the trivector part, without 
sources) d A [E /\ v/c + (IB) ■ v] =0. In contrast to all previous results, in 
the vector part, i.e., the part with sources, of the field equation there are two 
different velocities u and v. From these field equations it is concluded 3 that 
in the field equation with sources, (16) or (18) in Ref. 3, the usual Amper- 
Maxwell law and Gauss's law are inseparably connected in one law. Similarly, 
Faraday's law and the law that expresses the absence of magnetic charge are also 
inseparably connected in one law, the field equation without sources, Eq. (17). 3 
As shown in Sees. 6 and 7 in Ref. 3, this inseparability is an essential difference 
relative to the usual Maxwell equations with the 3-vectors. 

Next, we mention an important result regarding the usual formulation of 
electromagnetism, as in Ref. 9, which is presented in Ref. 10 and discussed in 
Refs. 11, 12 and 3. It is argued 10 that an individual vector has no dimension; 
the dimension is associated with the dimension of its domain. Hence, the time- 
dependent E(r,t), B(r,i), D(r, t) etc. cannot be the 3-vectors, since they are 
defined on the spacetime. Therefore, we use the term 'vector' for a geometric 
quantity, which is defined on the spacetime and which always has in some basis 



P = M-u/c, M = cI(MAu/c) 



(4) 



n =D Au/c+ (IH) - u/c 2 



(5) 
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of that spacetime, e.g., the standard basis {7^}, four components (some of them 
can be zero). Note that vectors are usually called the 4- vectors. However, an 
incorrect expression, the 3- vector, will still remain for the usual E(r,t), B(r,f), 
D(r, t) etc.. Moreover, recently, 13-17,11 it is proved that, contrary to the general 
belief, the UT of the 3-vectors of the electric and magnetic fields, E(r,t) and 
B(r,£) respectively, see, e.g., Eqs. (11.148) and (11.149) in Ref. 9, differ from 
the LT (boosts) of the corresponding 4D quantities that represent the electric 
and magnetic fields. As explained, 13 " 17,11 the fundamental difference between 
the UT and the LT of the electric and magnetic fields is that in the UT, e.g., the 
components of the transformed E' are expressed by the mixture of components 
of E and B, and similarly for B', Eq. (11.148). 9 The UT of the 3-vectors E and 
B are given, e.g., by Eqs. (11.149) 9 and they are 

E' = 7 (E + /3x c B)-(7 2 /(l + 7 ))/3(/3.E), 

B' = 7(B-(l/c)/3xE)-( 7 2 /(l + 7))/3(/3-B), (6) 

where E', E, /3 and B', B are all 3-vectors. All what is stated for the 3-vectors 
E and B and their UT holds in the same measure for the couple of the 3-vectors 
P and M and their UT 

P = 7 (P' + /3xM7 c )-(7 2 /(l + 7 ))/3(/3.p'), 
M = 7 (M'-/3xcP')-(7 2 /(l + 7))/3(/3-M'), (7) 

see the equations, e.g., Eq. (4.2), 18 or Eqs. (18-68) - (18-71), 19 or (6.78a) and 
(6.81a), 20 etc. 

However, the correct LT always transform the J^D algebraic object (vector, 
bivector) representing the electric field only to the electric field, and similarly 
for the magnetic field. 

In order to explain this fundamental difference between the LT and the 
UT let us introduce the frame of 'fiducial' observers as the frame in which 
the observers who measure fields E and B are at rest. That frame with the 
standard basis {7^} in it is called the 70-frame. In the 70-frame v — C70 and 
therefore E from (J2) becomes E — F ■ 70 and it transforms under the active 
LT (Eqs. (10) and (11) in Ref. 11) in such a manner that both F and the 
velocity of the observer v = C70 are transformed by the LT, see Eq. (12) 11 
(E = F ■ 70 — > E' = R(F ■ j )R = (RFR) ■ (Ry R)). As explained in Ref. 11, 
Minkowski, in Sec. 11. 6, 1 showed that both factors of the vector E, as the product 
of one bivector and one vector, has to be transformed by the LT. However, it 
is worth mentioning that Minkowski in all other parts of Ref. 1 dealt with the 
usual 3-vectors E, B, D, etc.. These correct LT give that 

£' = £ + 7 (£-/?){7o-(7/(l+7Mr, (8) 

Eq. (13). 11 In the same way vector B transforms and vectors P, M as well, but 
for P and M the LT, like ([5]). are the transformations from the rest frame of 
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the medium (u — C70). For boosts in the direction 71 one has to take in that 
equation that f3 — (3ji (on the l.h.s. is vector /3 and on the r.h.s. (3 is a scalar). 
Hence, in the standard basis and when /3 — f3-fi that equation becomes 

E w lv = -pjE 1 ^ + 7^ X 7i + E 2 l2 + E 3 l3 , (9) 

what is Eq. (14). 11 The most important result is that the electric field vector 
E transforms by the LT again to the electric field vector E' ; there is no mixing 
with the magnetic field B. The same happens with vectors P and M . 

On the other hand, if in the transformation of E = F ■ 70 only F is trans- 
formed by the LT R, but not the velocity of the observer v = c-fo (E — F--fo — !• 
E' F = (RFR) ■ 70, Eq. (15) 11 ), then, in the standard basis and when j3 — j3ji, 
one finds 

E%«f v = E x lx + 7 (£ 2 - cf3B 3 ) l2 + 7 (£ 3 + cpB 2 )-fy. (10) 

what is Eq. (17). 11 It is visible from the comparison of Eq. (fTU|) with Eq. 
(11.148) 9 that the transformations of components (taken in the standard basis) 
of E' F are exactly the same as the transformations of E XtVtZ from Eq. (11. 148), 9 
i.e., the components from ©. 



3. Constitutive Relations for Moving Media in Geometric Terms 

Let us consider the case of a simple medium with homogenous and isotropic 
electric and magnetic properties. In accordance with Sec. 2 we formulate the 
constitutive relations for moving media using the generalized magnetization- 
polarization bivector M. and the electromagnetic field bivector F 

M-u = EqXeF -u, (11) 

(IM) ■ u = { X b/^c 2 )u ■ (IF) (12) 

and the electric and magnetic susceptibility (xe, Xb)- 

Using the decompositions of F (p} and M (J3j we get from (fTT|) how the 
polarization vector P(x) depends on E, B and u, v, 

P = (e oXE /c){(l/c)[(u ■ v)E - (u ■ E)v] + {ufxvf\ B)I] (13) 

and from (|12[) how the magnetization vector M(x) depends on E, B and u, v, 

M = £oXs{[(w • v)B - (u ■ B)v] - (l/c)(u A v A E)I}. (14) 

Observe that P in (fT3"|) and M in ([T4^) contain both velocities u, the bulk velocity 
vector of the medium, and v, the velocity vector of the observers. The relations 
(fT3"|) and (fT4l are the basic relations that are obtained in this paper. In our 
approach, the relations (fT3)) and (|T4)) replace the constitutive relations with the 
3-vectors (|23|) and (|24|) . which are equivalent to Minkowski's equations (|22|) . 
The equations (|2"3")l . (|24l) and (|2^|) are given below in Sec. 5. 
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In the special case when v = u, an observer frame comoving with medium, 
[p]) and CE1} yield 

P = eoXEE, M = { X b/»o)B, (15) 

which are the familiar linear constitutive laws for the case of an electrically 
neutral, isotropic, non-dispersive, polarizable medium at rest. Introducing the 
relative permittivity e r , e r = 1 + XEj and the relative permeability fj, r , fi r = 
1/(1 — xb), the constitutive laws (fT5")) can be written as P = Eo(e r — 1)E, 
M = (l//io)(l — 1/ ' n r )B. In this case (v = u) the excitations D, D = eqE + P 
and H, H = (l//xo)-B — M can be introduced, as in Sec. 2, and the rest frame 
constitutive relations take the usual forms 

D = eE, H = (l//i)B, (16) 

where e = £o£ r and \i = /Jo/J r - 

The last term in ([TU)) and (fH)) describes the magnetoelectric effects in a mov- 
ing dielectric. According to the last term in (|13[) a moving dielectric becomes 
electrically polarized when placed in a magnetic field, the Wilsons' experiment. 21 
Similarly, according to the last term in (fT4)) a moving dielectric becomes mag- 
netized when it is placed in an electric field, Rontgen's experiment. 22 



4. Constitutive Relations for Moving Media in the Standard Basis 

The equations (fTTj) - (fit))) are all coordinate-free relations. In order to compare 
them with some usual formulations we have to represent all 4D geometric quan- 
tities from them in the standard basis {7^}. Then, in the {7^} basis, Eq. (fT3l) 
becomes 

P^m = (^Xe/c)[(1/c)(E^v v - E»v») + e^ af} v a B p \u vl ^ (17) 

whereas (fH)) takes the form 

M% = s oX b[(BV - B u v fi ) + (l/c)e^E a vp\u vllL . (18) 

We shall examine Eqs. (TT71) and (fTE)) taking that the laboratory frame, 
the S frame, is the 70-frame (v — C70, — (c, 0,0, 0), E° = B° = 0) in 
which the material medium, the S' frame, is moving with velocity u, u v — 
(luC, juU 1 , 'full 2 , 'f u U 3 ) , U k are the components of the 3-velocity U and j3 u = 
|U| /c. Then, in the laboratory frame, Eq. (fTT)) becomes 

P = e XElu[E l {U l /c) 1() + {E 4 + e^UjBM. (19) 

It can be seen from (JTHJ) that, e.g., if E» = (0,0,0,0), B^ = (0,0,0,-B 3 ), 
u ^ = {"fuCjjuU 1 ,0,0), the components of the polarization are 

P' i = (0,0,P 2 = e oX E7 u U 1 B 3 ,0) ] (20) 
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these components correspond to the 'translational' version of Wilsons' experi- 
ment. As stated on page 13, 23 the magnetoelectric effect 'in a moving dielectric 
is of course anisotropic in the sense that the polarization depends upon an ap- 
plied field which is perpendicular to the direction of motion and this polariza- 
tion is then perpendicular to both the applied field and the direction of motion.' 
Comparing P 2 from (f^Uj) with P y from Eq. (1.14) 23 (see also the term x"em) 
in Fig. 3.3 and Eq. (2.9) 23 ), we see that for the above conditions (the electric 
field is absent) these two expressions differ only in the j u factor; P 2 contains 
j u whereas P y contains j 2 . 

However, it is worth mentioning that the complete P y (both, the 3-vectors 
of the electric and magnetic fields exist) from Eq. (1.14), 23 or the components 
of the 3-vector P from Eq. (2.9) and Fig. 3. 3, 23 can be compared only with the 
spatial components of (fT9"|) , but P in (jT9"]) contains the temporal component as 
well. Again, the spatial components of (|19p and the complete P y differ only in 
the 7u factor. 

Similarly, from (|18[) , we find 

M = ( X B/^hu[B' l (Uyc) l0 + (B l - e^UjEk/c 2 )^}. (21) 

In the case of low velocities of the medium, f3 u <C 1, the equations flT9"|) and 
(|2"Tj) are almost unchanged to first order of f3 u ; only it is taken that j u ~ 1. 
Observe that in flT9"|) and (|2"Tj) the terms with 70 are comparable to the terms 
containing UjBkji and (l/c 2 )UjEkji, respectively. Thus, to first order of /3 U , 
the terms P° and M° cannot be neglected relative to the complete P % and M l , 
respectively. 

The constitutive relations (fT9|) and (|2"Tj) significantly differ from all previous 
constitutive equations for such simple moving media. They will be compared 
with the usual formulations with the 3-vectors. 



5. Comparison with Minkowski's Constitutive Relations 
and with Their Derivations 

Let us examine some of the usual derivations of the constitutive relations for 
moving media. As already mentioned, almost the entire physics literature deals 
with the electromagnetic excitation tensor H, i.e. with the electric D and mag- 
netic H excitations and the constitutive relations refer to the connections be- 
tween them and F, i.e. E and B, respectively. 

Minkowski 1 was the first who derived the constitutive relations for the 3- 
vectors 

D + (l/c 2 )UxH = e(E + UxB) 

B-(l/c 2 )UxE = /i(H-UxD), (22) 

Eqs. (C) and (D) in Sec. 8. 1 From that time, the equations ([22]) are considered 
to be the fundamental constitutive equations for moving media in the whole 
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physics community. They are subsequently derived in different ways and used 
in numerous papers and textbooks on the electromagnetism. 

5.1 Comparison with Minkowski's constitutive relations that are 
obtained by means of the 3-vectors 24 and 
in the covariant approach 25 

One way of the derivation of Eqs. (|22l) is presented in Ref. 24. There, these 
equations, i.e. Eqs. (7), 24 are derived using Minkowski's crucial hypothesis that 
the relations with the 3-vectors D = eE, H = (l//i)B, which correspond to (|TB)) . 
retain their form in the moving frame S' with the same e and //, i. e., D' = eE', 
H' = (l//i)B'. Then, in Ref. 24, the UT © and the similar UT for D' 
and H' (Eqs. (3) 24 ) are used to obtain Minkowski's constitutive equations (|22[) . 
Furthermore, in Ref. 24, the constitutive relations, which explicitly express D 
and H in terms of E and B, are derived (Eqs. (10) 24 ) from (|22p . Introducing 
D = e E + P and H = (l//Xo)B - M into Eqs. (10) 24 we find the constitutive 
equations which explicitly express P as a function of E and B 

P - 7 2 £o{xb[E-c- 2 U(U-E)+UxB] 

+Xb[Ux (B-c- 2 UxE)]} (23) 

and also M in terms of E and B 

M = ( 7 7w>){xb[B- C - 2 U(U-B)- c - 2 UxE] 

-c- 2 Xe {Ux (E + UxB)]}. (24) 

The equations (|2"3")l and (f2~4"|) are equivalent to Minkowski's equations (|2"2"|) . 

Introducing the relative permittivity e r and the relative permeability /i r 
instead of the susceptibilities xe and \b the equations (f23|) and (|24|) become 

P = eo(s r — 1)E + -f 2 eo(e r — l//i r )U x (B — c~ 2 U x E), 
M = - l//i r )B -7 2 e (e r - l//i r )U x (E + U x B). (25) 

and they are, in the same way as Eqs. (f23)) and (f24)). equivalent to Minkowski's 
equations (|22|) . 

Similarly, in Ref. 25, the relations (J22|) (Eqs. (76. 9) 25 ) are obtained by the 
covariant generalization (only components, implicitly taken in the standard ba- 
sis) of the relations with the 3-vectors D = eE, H = (l//x)B, which correspond 
to our Eq. ((T6)) . The covariant generalizations from Ref. 25 are T-L x ^u^ L = 
eF^u^ ((76. 7) 25 ) and *F x ^ Ufl = ^U^u^ ((76.8) 25 ), where *F X ^ and *n Xf * 
are the dual tensors, i.e., only the components (*F a P = (l/2)e a/37S F 1 s)- Then, 
Minkowski's constitutive equations (f22|) with the 3-vectors are derived by Minkowski's 
identifications, 1 of components of F a P with components of the 3-vectors E 
and B {E l = F l ° , B l = (l/2c)e zjk0 F jk , see, e.g., Eq. (11.137) 9 ) and by 
the similar identifications for H al3 and the 3-vectors D and H {D l = H lQ , 
W = (c/2)e^ k0 n jk ). 
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However, as discussed, e.g. in Sec. 2, 11 the components are not the whole 
physical quantity. The mentioned identifications of components are synchroniza- 
tion dependent and the UT ^ (and the similar ones for D' and H') significantly 
differ from the LT (JgJ. 

On the other hand, our equations (fTTj) and (TT2"]) deal with coordinate-free, 
4D geometric quantities. Instead of Minkowski's identifications of components, 
and the same ones in Ref. 25, the mathematically correct decompositions of F 
|T]) and M ([3]) are used for the derivation of the constitutive equations ([TB"]) and 
(and, in the standard basis, (|TT)) and ([T5]) ) from (fTTl) and (fT2"]) . respectively. 

As mentioned above, the equations (|2"3"]) and (IM1) . or (j2"5)) , which are equiv- 
alent to (|2"2"|) , can be compared with Eqs. (fH?)) and (j2"Tj) . There are important 
differences between them. 

1) The equations (fT9|) and (|2"Tj) are with correctly defined 4D quantities that 
properly transform under the LT, whereas it is not the case with Eqs. (|23p and 
(12"4")) . or (|23|) . with the 3- vectors that transform according to the UT. 

2) Furthermore, (fH)]) and (|2"Tj) contain the term with 70, which cannot exist 
in the approach with the 3- vectors, i.e. in (|2"5|) and (|2"4"|) . or (|25p . 

3) Also, in (jT9j) , the polarization vector P contains only the electric suscep- 
tibility xe and not xb, whereas, in (|23l) . the polarization 3-vector P contains 
both susceptibilities xe and xb, (or, in ([2"5]). both, e r and Similarly, in (l2"Tj) , 
there is only x# and not X£, whereas, in l[2"4"]). there are both susceptibilities 
Xb and (or, in (|25j) . both, /x r and e r ). It can be seen from the derivation 
of (f23|) and (|24|) . or (|25[) . that both susceptibilities appear in the expressions 
for P ((231) and M ([24]) because the UT © and the similar ones for D' and H' 
(Eqs. (3) 24 ), or, equivalently, the UT ((TJ), are used and in them, e.g., in ((6]), E' 
is expressed by the mixture of E and B, and similarly, in ((7|), P' is expressed by 
the mixture of P and M. 

5.2. Comparison of the low velocity limits 

These differences remain in the low velocity limit, j3 u -C 1, as well. That limit 
is already examined for vectors P and M from (|19[) and (|2"Tj) . respectively. Only 
7„ ~ 1 is taken in (Jl9]) and ([21]) . 

In Ref. 24, the low velocity limit, i.e., the quasi-static approximation, is ob- 
tained in the same way by putting j u ~ 1 in Eqs. (10) for D and H expressed 
in terms of E and B, which led to Eqs. (II). 24 They, Eqs. (11), 24 are com- 
monly used in literature. In the formulation with the 3-vectors P and M the 
equations (|25|) are derived from Eqs. (10). 24 Therefore, in that approximation, 
the usual 3-vectors P and M are given by ([23]) but with j u ~ 1. This means 
that again both e r and /i r , i.e., both susceptibilities xe and Xb, appear in the 
quasi-static approximation in the expressions for P and M. In Ref. 24, it is 
argued that such a procedure with 7„ ~ 1 does not give a correct quasi-static 
approximation, because the obtained equations, (ll), 24 'do not obey the group 
additivity property' Hence, according to Ref. 24, the equations (ll) 24 have 
to be replaced by two well-defined Galilean limits of Minkowski's constitutive 
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equations f(|22| here), the magnetic and electric limits, i.e. with two sets of 
low-velocity formulae. These two sets are Eqs. (12) and (13) in Ref.24 and it is 
stated there: 'Both sets (12) and (13) do form a group.' 

However, the UT © and © and the similar UT for D' and H' (Eqs. (3) 24 ) 
are not the LT and therefore constitutive equations with the 3- vectors, (|2"2"j) . 
or Eqs. (10), 24 , and also (|2"3"|) and (j2~4")l . or (|2l)|) . here, are not the relativistic 
constitutive equations and also Eqs. (12) and (13) 24 are not a quasi-static 
approximation of the relativistically correct constitutive relations. Thus, as 
stated above, all differences from Sec. 5.1 remain in the low velocity limit as 
well. 

In the special case of moving (j3 u <C 1) nonmagnetic (fi r = 1, i.e. \b = 0) 
media the 3D polarization P from ([25]) . to first order of j3 u , becomes 

P=£oXb(E + UxB). (26) 

The equation (gBJ) is Eq. (9-19) (2), 19 or the last equation in Problem 6.8. 20 

But, even for (3 U <C 1, P from (|19p will have the same form for both cases 
Xb / and xb = 0, because it does not depend on xb (P — £qXe[E 1 {U 1 / c)lo + 
(E i + £ Oljk U B k )^i]). As already stated, to first order of (3 U , P° cannot be 
neglected relative to P l , which means that we have not an equation that would 
correspond to |(25ty . 

5.3. Comparison with constitutive equations that are obtained 
by means of exterior forms 2 

It is worth noting that Minkowski's constitutive equations (|22|) are also obtained 
in Ref. 2 using exterior calculus, i.e., with the use of abstract 4D geometric quan- 
tities, Eqs. (E.4.28) and (E.4.29). 2 At that place (p. 353) Hehl and Obukhov 
stated: 'Originally, the constitutive relations (E.4.28), (E.4.29) were derived by 
Minkowski with the help of Lorentz transformations for the case of a flat space 
time and a uniformly moving medium. We stress, however, that the Lorentz 
group never entered the scene in our derivation above.' It can be concluded 
from these statements that the authors 2 believe, as all others, that the UT ([6]) 
are the relativistically correct LT. However, as discussed in Sec. 2, the LT are 
given by Eq. (jHJ and not by ([6]). Furthermore, in Ref. 2, the expressions for the 
polarization P and the magnetization M, the equations (E.4.30) and (E.4.31), 
respectively, are derived from the constitutive relations (E.4.28), (E.4.29). The 
expressions for P and M, (E.4.30) and (E.4.31), respectively, from Ref. 2, 
are very similar to the equations with the 3-vectors, (|23|) and (|24[) . respectively, 
which are obtained from (f^^j) . Indeed, all quantities in the constitutive relations 
(E.4.28), (E.4.29) and in the expressions for P and M, (E.4.30) and (E.4.31), 
respectively, from Ref. 2, are in the 3D subspace, i.e., like the usual 3-vectors. 
(Observe that v that enters into these relations is the 3-velocity 1-form, (E.4.27), 
or the 3-velocity vector. Also, the Hodge star operators in these expressions are 
defined by the 3-space metrics of the laboratory and material foliations.) This 
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means that these results 2 also significantly differ from the constitutive relations 
(fT3|) and (fT4|) in which all quantities are the 4D geometric quantities and there 
is no the space-time split. Let us describe in more detail the calculation from 
Ref. 2, which led to the mentioned results. 2 

In Ref. 2, Minkowski's constitutive relations (E.4.28), (E.4.29) and the ex- 
pressions for P and M, (E.4.30) and (E.4.31), respectively, are derived using 
(1+3) - splitting of spacetime both in the laboratory frame and in the frame of 
moving macroscopic matter, the material frame. The starting relations in that 
derivation are the decompositions ((1+3) - splitting) of the electromagnetic 
excitation tensor T-L and of F in the laboratory frame (unprimed quantities), 
(H = - H A dt + D, F = E A dt + B) Eq. (E.4.14), and in the material frame 
(primed quantities), (H = - H' A dt' + D' , F = E' A dt' + B') Eq. (E.4.15). 
(In both equations our notations is used.) There is an assumption in connec- 
tion with the decompositions (E.4.14) and (E.4.15) in Ref. 2. It is: 'Clearly, 
we preserve the same symbols T-L (our notation) and F on the left-hand sides 
of (E.4.14) and (E.4.15) because these are just the same physical objects. In 
contrast, the right-hand sides are of course different, hence we use primes.' 

From the mathematical point of view there are no reasons for such argu- 
mentations. If H , dt, D, E, B arc transformed then T-L and F have to be 
transformed as well, as can be concluded going to a uniformly moving medium, 
i.e., when the LT are applicable. Let us explain these assertions in more detail. 
There is no mathematical procedure by which in one equation some parts of it 
are transformed and the other parts are not. It is not possible that there is one 
law for the transformations of some J^D geometric quantities and another law 
for the transformations of the other quantities. As mentioned, e.g., in Ref. 11, 
"any multivector M transforms by the active LT in the same way, Eq. (II), 11 
M — > M' = RMR, where the Lorentz boost R is given by Eq. (10). 11 Hence, if 
R is applied to (p}, R(F = E A v/c + (IcB) -v/c)R, then F is also transformed, 
together with E, B and v. Thus, the assertion that F is unchanged under the 
active LT is without any mathematical justification. 

In addition, let us compare the relations (E.4.14) and (E.4.15) with the 
consideration from Ref. 26. In Sec. 4 'Premetric electrodynamics in exterior 
calculus' in Ref. 26, it is argued that the decomposition of H, Eq. (31), 26 if 
put in matrix form, yields the identification (6)i in Ref. 26. Furthermore, the 
field strength F is decomposed in the same way as T-L, i.e. 'into two pieces: 
one along the 1-dimensional time t and another one embedded in 3-dimensional 
space (a, 6=1,2, 3.).' Thus, in the notation, 26 F = EAdt + B = E a dx a Adt + 
(l/2)B a bdx a Adx b , Eq. (34). 26 If put in matrix form, it yields the usual identifi- 
cation (6) 2 in Ref. 26, i.e., in our notation, E* = F l °, B l = (l/2c)e ijk0 F jk . The 
decompositions (31) and (34) 26 are the same as the decompositions (E.4.14). 2 
As already stated, in the material frame, the same decompositions of H and 
F (F — E' A dt' + B') are performed, (E.4.15), 2 and they yield the usual 
identifications (6)1 and (6)2 in Ref. 26, but now with primes, E n — F n0 , 
B n — (l/2c)e^ kQ Fj k . This discussion reveals that the (1+3) - splitting of space- 
time, i.e., the usual identifications (6)1, (6)2 in both frames, are the real cause 
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that, although the authors 2 worked with exterior forms in the 4D spacetime, 
they obtained the same results for Minkowski's constitutive relations and for D 
and H ((E.4.25) and (E.4.26) 2 ), i.e., for P and M, as those that were found by 
means of the 3- vectors and their UT, e.g., in Ref. 24. 

However, as already discussed, e.g., in Refs. 27 and 11, and also in Ref. 12, 
the (1+3) - splitting of spacetime and the usual identifications of components of, 
e.g., F a/3 (implicitly taken in the standard basis) with components of the 3D E 
and B in both frames are synchronization dependent and they are meaningless 
in the V synchronization. We note that different synchronizations are nothing 
else than different conventions and physics must not depend on conventions. (In 
Ref. 26, Hehl (Eq. (6)2) quotes Minkowski's identifications from Sec. 3 1 and 
states: 'We stress that the identifications in (6) are premetric, they are valid on 
any (well-behaved) differential manifold that can be split locally into time and 
space.' This means that Hehl considers that 1+3 decomposition, i.e., the space- 
time split, is also - premetric. But, if a nonstandard basis, the {r^} basis with 
the 'r' synchronization is used, i.e., when the appropriate metric is used, then 
it is not possible to make the usual identifications, that is, 1+3 decomposition, 
i.e., the space-time split. That metric is discussed in, e.g., the text below Eq. 
(14) 12 : '.. the components g^, r of the metric tensor g a j, are ga^ r — 0, and all 
other components are = 1.' Hence, both, the usual identifications (Eqs. (6) 26 ), 
and 1+3 decomposition, i.e., the space-time split, are meaningful only when the 
Minkowski metric, e.g., diag(l, —1, —1, —1), is used. Thus, these identifications 
depend on the chosen metric and therefore they are not premetric. 

(For the {r M } basis see, for example, Refs. 27, 12.) In the {r^} basis, 
it holds that F r 10 = E 1 + cB 3 - cB 2 . The relation for F r 10 shows that the 
components have no definite physical meaning since they are dependent on the 
chosen synchronization. Only in the {7^} basis does it hold that E l = F l °, 
p t _ M io^ etc _ According to that, the usual 3-vectors E, B, D, H, P, M, 
etc., where, e.g., P =M 10 i + M 20 j + _M 30 k (i, j, k are the unit 3-vectors) have 
no definite physical meaning, since the components A4 M are dependent on the 
chosen synchronization.) 

Comparing, for example, the decompositions of F from Ref. 2 and 26 (the 
above quoted Eq. (34) 26 ) with our decomposition F = E A v/c + (IcB) ■ 
v/c, Eq. ([T]), i.e., in some arbitrary basis {e^} (it can be the standard basis 
{7^}, or the {r M } basis with the 'r' synchronization, etc.) F = (5 a/S '^E^v" + 
c£ a l 3 ^ Ve ^B e l/ )e a Aej3 (a, f3,fj,,v = 0, 1, 2, 3 and E£, B^, v£ are the components 
in the {e^} basis), we see that there are important differences between them. 
All quantities, including v, in our relations are the 4D quantities and there is no 
(1+3) - splitting of spacetime, which means that they are more general than the 
decompositions from Ref. 2 and 26. If, for example, the magnetization 1-form 
M, Eq. (E.4.31), 2 would be written in the {7^} basis then it would contain 
in the material frame and in the laboratory frame only spatial components and 
not the time component in contrast to the relation (|2"TT) (and (|T5)) V Observe 
also that relations (JT2J) and (IT4")) . or in the {7^} basis (JTTJ) and (fT5)l . contain 
both the velocity of the observer v and the velocity of the considered medium 
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u, whereas there is only one velocity, the velocity of the considered medium, in 
the equations (E.4.25)-(E.4.31). 2 

Furthermore, the obvious similarity of the results with exterior forms 2 and 
the usual results with the 3- vectors and their UT, e.g. from Ref. 18 (or 24), can 
be nicely seen comparing Sec. E.4.5 under the title 'The experiments of Rontgen 
and Wilson & Wilson' in Ref. 2 and Sees. 5.2 and 5.3 in Ref. 18. It is visible 
from Eq. (E.4.49) and Fig. E.4.2 2 that in the matter- free region, region (2), 
there is only an electric field E(o), whereas there are both the electric field Ern 
and the induced magnetic field B/u in the slab moving with constant velocity. 
Completely the same result for that case is obtained in the third equation in 
(5.18) in Sec. 5. 2, 18 where it is exclusively dealt with the 3-vectors and their 
UT. The same conclusion holds for the comparison of the treatments of the 
Wilsons' experiment, Sec. E.4.5 2 and Sec. 5.3. 18 

From the above discussion it is visible that either by using the 3-vectors and 
their UT, or by using the modern mathematical language, the exterior forms, 
always, in all previous treatments, it is obtained that, e.g. the electric field in 
one frame is 'seen' as slightly changed electric field and the induced magnetic 
field in the frame relatively moving with constant velocity. The same holds for 
the polarization and the magnetization and also for the electric and magnetic 
excitations. 



6. The Physical Explanation of the Magnetoelectric Effect 
in Moving Media 

The consideration in the previous sections shows that there are important dif- 
ferences in the form of Minkowski's constitutive relations (|2"2"j) . or equivalently 
the constitutive relations for P (f23|) and M (l24l) . i.e., (|25l) . and the constitutive 
relations for vectors P (Tl3l) and M (fl4|) . However, there is also a significant 
difference in physical interpretation of these constitutive relations. Particularly, 
this refers to the interpretation of the magnetoelectric effect for moving media. 

In previous approaches it is considered that the term which describes how the 
magnetic field influences the polarization {^ 2 e a (e r — l//i r )U x B in Eq. (|2"5"j) ) 
and the term which describes how the electric field influences the magnetization 
(-7 2 e (£r - 1/V)U x E in Eq. (J25J) ) are determined by the UT © for E and 
B. Namely, according to the UT (|6]), in the rest frame of a moving medium the 
magnetic field 3- vector from the laboratory frame is 'seen' as a slightly changed 
magnetic field 3-vector and an induced electric field 3-vector. That induced 
electric field 3-vector interacts with the electric dipole moments 3-vectors by 
the interaction term E • d giving the polarization P as a 3-vector. Similarly, 
in the rest frame of a moving medium the electric field from the laboratory 
frame is 'seen' as a slightly changed electric field and an induced magnetic field. 
That induced magnetic field interacts with the magnetic dipole moments by 
the interaction term B • m giving the magnetization M (all are 3-vectors). In 
my opinion, it is very strange that, e.g., the magnetic field from a permanent 
magnet is 'seen' as an induced electric field in a relatively moving frame. What 
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is the physical reason which causes that one field becomes another field in a 
relatively moving frame? How can it be that a mere uniform moving transforms 
one field to another one? It cannot be that the transformations describe the 
physics, but, on the contrary, the physics has to describe how to get the correct 
transformations of the fields. 

If the electric and magnetic fields are interpreted as vectors defined on the 4D 
spacetime then it is very natural to have that a magnetic field vector remains 
the magnetic field vector in a relatively moving frame. There is no physical 
cause that can change one field to another one in a relatively moving frame. 
According to the LT ©, an electric field vector transforms again to the electric 
field vector and similarly for a magnetic field vector. Hence, the explanation 
for the magnetoelectric effect in moving media is completely different in our 
approach. 

Instead of dealing with the electric and magnetic dipole moments 3-vectors 
d and m we deal with the dipole moment bivector D, as a primary quantity for 
dipole moments, which does have a definite physical reality. The decomposition 
of D into the electric and magnetic dipole moments vectors d and m, respec- 
tively, and the unit time-like vector u/c, and the expressions for d and m, which 
are obtained from D and determined relative to u are 

D = d A u/c + (ml) ■ u/c 2 , 

d = D-u/c, m = cI(D A u/c); (27) 

compare with © and ((!]). In this case the vector u is the bulk velocity vector 
of the medium as in ([3]) and Q • The interaction term in the Lagrangian for the 
interaction between the electromagnetic field F and the dipole moment bivector 
D can be written as a sum of two terms 

L mt = F-D = (l/c 2 )[-(E-d + B-m)(v-u) + (E-u)(vd) (28) 
+(B ■ u)(v ■ m)] + (l/c 3 )[(E Am- c 2 B Ad) Av A u]I. 

In the tensor formulation, the relations (f2~T|) and (|2"51) are given by Eqs. (2) 
and (3), respectively, in Rcf. 12 (they are first reported in Ref. 28). Observe 
that every term on the r.h.s. of (|28| contains both velocities u and v. As 
seen from the last two terms they contain the direct interaction of E with m, 
and B with d. These terms give the physical explanation for the existence of 
the magnetoelectric effect in moving media. Moreover, there is no need for any 
transformation. We only need to represent E, d, B, m, u and v from (1281) in 
the standard basis and then to choose the laboratory frame as our 70-frame. It 
can be seen from the discussion of Eq. (25) 12 that in the laboratory frame, as 
the 70-frame, and in the low velocity limit, we can neglect the contributions to 
Li n t from the terms with dP and m°; they are u /c 2 of the usual terms E ■ d or 
B ■ to. Then, what remains from (1281) is 

L^t = -((£«**) + [B im 1 )) - (l/c 2 )e°^ k (E imk ~ c 2 B. l d k )u J . (29) 

This is, to order 0(u 2 /c 2 ), relativistically correct expression with vectors for Li nt . 
The last two terms that contain the direct interactions between E and to and 
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between B and d are not taken into account in any of the previous investigations 
of the magnetoelectric effect in moving media. They are u/c of the usual terms 
with the direct interaction of E with d and B with m. (The expression (|2"9"|) is 
first reported in Ref. 29 in connection with the EDM searches.) 



7. A Brief Discussion of the General Constitutive Relations 

In this paper, the consideration is restricted to the constitutive relations and the 
magnetoelectric effect in moving media with homogeneous and isotropic electric 
and magnetic properties. For them, the rest frame constitutive relations are 
given by Eqs. (|16|) . The general constitutive relations which cover dielectric, 
magnetic and magnetoelectric behavior are considered, e.g., in Ref. 26 and in 
more detail in Refs. 30, 2, 31 and 23. Mainly, e.g., Refs. 2, 26, 30, 31, these 
more general constitutive relations link the electromagnetic excitation H and F, 
Ha/3 — K Q( g 7<5 F 7 5 (our notation), Eq. (17), 30 where l K a/ 3 lS (x) is the twisted con- 

2 1 

stitutive tensor of type ' or, in another form, e.g., *U afi = x aM& F 1& ( our 

notation), Eq. (21), 26 where *U afi = (l/2)e a ^ 5 H-ys, X at>1& is 'a constitutive 
tensor density of rank 4 and weight +1, with the dimension [x]=l/resistance.' 
In Refs. 23, the general constitutive relations are established by expressing Ai 
as a linear function of the electromagnetic field F, nocMafs = (l/2)£2g-f7<5 (our 
notation), Eq. (2. 45), 23 where '£^1 is the general susceptibility tensor, which 
is dimensionless because of the choice of the constant /xoc.' However, in all 
these treatments, all quantities in the equations are not geometric quantities 
but components in the standard basis, which means that the components of F, 
H, M are obtained by the usual identifications, i.e., they are considered to be 
the components of the 3-vectors E and B, D and H, P and M, respectively. 
This can be nicely seen from Eqs. (2.47) and (2.48) 23 and from their comparison 
with Eq. (2. 9). 23 Hence, these relations are not so general relations and they 
are not premetric, as stated, e.g., in Ref. 26, because, as already mentioned in 
Sec. 5.3, the usual identifications and the space-time split are meaningless, e.g., 
in the basis with the 'r' synchronization. 

In contrast with the usual covariant approach with coordinate-dependent 
quantities, all relations flTJ - ©, ©, dTTJ) - (JlgJ) an d (27]) - (28| are coordinate- 
free relations, which are written in terms of the abstract 4D geometric quantities. 



8. Conclusions 

The constitutive relations for P (fT3"|) and M (fT4l contain both velocity vectors 
u and v and thus they differ from all previous expressions. They arc formulated 
in terms of coordinate-free quantities that correctly transform under the LT 
(JSJ, whereas, as explained in Sees. 5 - 5.2, it is not the case with Minkowski's 
constitutive relations ([23]) and (|24|) . or (|25[) . with the 3-vectors that transform 
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according to the UT © and ([7]), which are not the LT. As discussed in Sec. 
5.3, Minkowski's constitutive relations for P and M, (E.4.30) and (E.4.31), 
respectively, from Ref. 2, are obtained by the (1+3) - splitting of spacetime and 
consequently they are equivalent to the usual expressions with the 3-vectors 
(|2"3")1 and Hence, these relations, 2 which are obtained using exterior forms, 
also differ from our results (jT3"]l and (|T4")) . The differences that are quoted 
at the end of Sec. 5.1, points 1) - 3), could be used for the experimental 
examination and comparison of the results presented here and the constitutive 
relations which are obtained in the usual formulations either with the 3-vectors 
or with exterior forms. Furthermore, a completely new physical explanation 
of the magnctoelcctric effect in moving media that is presented in Sec. 6, Eq. 
|28|) , or Eq. (f29| , offers the possibility for the experimental investigations of the 
magnetoelectric effect from the relativistically correct point of view. Regarding 
the importance of the magnetoelectric effect the results obtained in this paper 
could be enough important in different physical applications. 
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